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ABSTRACT 

I study the evolution of haloes density profiles as a function of time in the SCDM 
and ACDM cosmologies. Following Del Popolo et al. (2000) (hereafter DP2000), I 
calculate the concentration parameter c = r„/a and study its time evolution. For a 
given halo mass, I find that c(z) oc l/{l + z) both in the ACDM and SCDM cosmology, 
in agreement with the analytic model of Bullock et al. (1999) (hereafter B99) and N- 
body simulations. In both models, a(z) is roughly constant. The present model predicts 
a stronger evolution of c{z) with respect to Navarro, Frenk & White (1997) (hereafter 
NFW97) model. Finally I show some consequences of the results on galaxy modelling. 
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1 INTRODUCTION 

The structure of dark matter haloes is of fundamental importance in the study of the formation and evolution of galaxies and 
clusters of galaxies. From the theoretical point of view, the structure of dark matter haloes can be studied both analytically 
and numerically. A great part of the analytical work done so far is based on the secondary infall model (SIM) introduced by 
Gunn & Gott (1972). Calculations based on this model predict that the density profile of the virialized halo should scale as 
p oc r^^^'^. Self-similar solutions were found by Fillmore & Goldreich (1984) and Bertschinger (1985), who obtained a profile 
of p (X r~^'^^ . Hoffman & Shaham (1985) (hereafter HS) considered a scale-free initial perturbation spectra, P{k) oc fc". They 
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2 A. Del Popolo 

showed that p oc r~" with a — ^^^p^, thus recovering Bertschinger's (1985) profile for n = and f2 = 1. They also showed 
that, in an open Universe, the slopes of the density profiles steepen with increasing values of n and with decreasing SI, reaching 
a profile p cx r~* for f2 0. 

N-body simulations, such as that of Quirm, Salmon & Zurek (1986), West, Dekel & Oemler (1987) and Efstathiou et al. 
(1988) arrived at conflicting results implying that better numerical resolution were needed to settle the issue. Recent results 
from higher resolution simulations (Navarro, Frenk & White 1995, 1996, 1997 (hereafter NFW95, NFW96, NFW97); Lemson 
1995; Cole & Lacey 1996; Tormen, Diaferio & Syer 1998) obtained, using different codes and different setups for the initial 
conditions, agreement in the conclusion that halo density profiles do not follow a power law but develope a universal profile, 
a one parameter functional form that provides a good fit to haloes over a large range of masses and for any scenario in which 
structures form due to hierarchical clustering, characterized by a slope (3 — — ~1 near the halo centre and /? = —3 at 

large radii. In that approach, density profiles can be fitted with the functional form: 
p('')nfw _ 1^11 

a \ a } 

where pb is the background density and 5n is the central overdensity [below I shall refer to equation (|l|) (NFW97) as the NFW 
profile]. The scale radius a, which defines the scale where the profile shape changes from slope /3 < —2 to /3 > —2, and the 
characteristic overdensity, 5n, are related because the mean overdensity enclosed within the virial radius is ~ 180. I recall 
that according to NFW96, 97), a is linked to a "concentration" parameter, c, by the relation a = ^ and the parameter c is 
linked to the characteristic density, Jn, by the relation: 
_ 200 

3 ln(l + c)-c/(l + c) ^ ' 

The scale radius and the central overdensity are directly related to the formation time of a given halo (NFW97). The power 
spectrum and the cosmological parameters only enter to determine the typical formation epoch of a halo of a given mass, and 
thereby the dependence of the characteristic radius, a, or the overdensity (5n on the total mass of the halo. (5n increases for 
decreasing virial mass, A/v. A natural reason for the fact that low-density haloes tend to show higher densities is that they 
typically collapse earlier, when the universe was denser. To model this trend, NFW97 proposed a step-by-step calculation of the 
density profile assuming that the characteristic density, is proportional to the density of the universe at the corresponding 
collapse redshift, z^. This model successfully predicts the (5n — Mv relation for different cosmological models at z = 0. The 
model has been also extended in NFW97 to predict the redshift dependence of the halo profile parameters but, as shown by 
B99 for the ACDM cosmology, the evolution of c (and consequently of 5n) is much stronger than in the NFW97 model. 
In this paper, I use the improved SIM introduced by DP2000 to determine c(^z) for both SCDM and ACDM models and to 
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Profiles of dark-matter haloes at high redshift 3 

compare it with prediction of NFW97 and B99 models. 

The plan of the paper is the following: in Section 2, I introduce the model. In Section 3, I show the results of the model and 
Section 4 is devoted to the conclusions. 

2 TIME EVOLUTION OF THE CONCENTRATION PARAMETER 

The simplest version of SIM considers an initial point mass, which acts as a nonlinear seed, surrounded by a homogeneous 
uniformly expanding universe. Matter around the seed slows down due to its gravitational attraction, and eventually falls 
back in concentric spherical shells with pure radial motions. The assumptions of SIM that are most often questioned are 
the spherical symmetry and the absence of peculiar velocities (non-radial motions): in the "real" collapse, accretion does not 
happen in spherical shells but by aggregation of subclumps of matter which have already collapsed; a large fraction of observed 
clusters of galaxies exhibit significant substructure (Kriessler et al. 1995). Motions are not purely radial, especially when the 
perturbation detaches from the general expansion. Nevertheless the SIM gives good results in describing the formation of dark 
matter haloes, because in energy space the collapse is ordered and gentle, differently from what is seen in N-body simulations 
(Zaroubi, Naim & Hoffman 1996). As I showed in a recent paper, DP2000, the discrepancies between the SIM and some 
high resolution N-body simulations are not due to the spherical symmetry assumption of the SIM but arises because of some 
non-accurate assumptions used in its implementation. As I showed in DP2000, the predictive power of the SIM is greatly 
improved when some problems of the previous implementations are removed. 

To begin with, the conclusion p oc for n < —1, claimed by HS, is not a direct consequence of the HS model, but it is 
an assumption made by the quoted authors, following the study of self-similar gravitational collapse by Fillmore & Goldreich 
(1984). In fact, as reported by the same authors, in deriving the relation between the density at maximum expansion and the 
final one, HS assumed that each new shell that collapses can be considered as a small perturbation to the gravitational field 
of the collapsed halo. This assumption breaks down for n < —1. 

Secondly, the assumption made by Hoffman & Shaham (1985) that S(r) oc ^(r) oc r-(='+"' is not good for regions internal 
to the virial radius, (see Peebles 1974; Peebles & Groth 1976; Davis & Peebles 1977; Bonometto & Lucchin 1978; Peebles 

3(3 + r.) 

1980; Fry 1984). In the inner regions of the halo, scaling arguments plus the stability assumption tell us that ^(r) oc r (f>+") , 
and we expect a slope different from that of HS. In other words, HS's (1985) solution applies only to the outer regions of 
collapsed haloes, and consequently the conclusion, obtained from that model, that dark matter haloes density profiles can be 
approximated by power-laws on their overall radius range is not correct. It is then necessary to introduce a model that can 
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make predictions also on the inner parts of haloes. 

Thirdly, according to Bardeen at al. (1986), (hereafter BBKS), the mean peak profile depends on a sum involving the initial 
correlation function, ^(r) oc r"'^+"', and its Laplacian, V^C('") « r~'^+"' (BBKS; Ryden & Gunn 1987): 



7'e(0 + ^V^e(r) 



C(0)V2 ^(1-^2) 
where v is the height of a density peak 

a{R,z) 

The variance a-{R, z) is given by: 



a^iR,z) ^ D'\z,Q.)^ J dkk'^P{k)W^{kR) (5) 

where the function D{z,fl) describes the growth of density fluctuations (Peebles 1980) and W(kR) is a top-hat smoothing 
function: 

W{kR) = (sin kR ~ kR cos kR) (6) 

7 and R, are two spectral parameters given respectively by: 
f k'^P(k)dk 

7 = — (7) 

[/ k'^P{k)dk J k^P{k)dk\ ' 



m-"^ (3) 



(4) 



R. 



3 f k^P(k)dkY'^ 

^ ^ ' (8) 



/ k^P{k)dk 

while j?(7J^, 7) is: 

3(1 - 7^) + (I.2I6 - 0.97") exp 

^(jy7,7) = _ 



(i)( 



3(l--72) + 0.45+(f^)' 
I recall that the z dependence of S is: 



1/2 

^ 2 



(9) 



s{z) = SoD(z,n) (10) 

being So the overdensity as measured at current epoch to. As can be seen for example in the case of a scale- free density 
perturbation spectrum (DP2000, equation (20)), the initial mean density obtained using the model of this paper is extremely 
different from that obtained and used in IfS. 

The flrst step to get c{z) is to calculate 5{r) for a given cosmology starting from the related spectrum. In order to calculate 
5{r) in the SCDM cosmology {fl = 1, h — 0.5, n = 1), I use the spectrum given by BBKS: 

P{k) = Ak'^ [In (1 + 4.164fc)]^ 

X (192.9 + 1340k + 1.599 x 10^ A:^ + 1.78 x 10^ A:^ + 3.995 x W^k*)~^^^ (11) 
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Profiles of dark-matter haloes at high redshift 5 

normalized by imposing that the mass variance at 8^~^Mpc is as = 0.63. For the ACDM model {Qni = 0.3, SIa — 0.7, h — 0.7), 
I also use the BBKS spectrum normalized as erg = 1. Supposing that energy is conserved, the shape of the density profile at 
maximum of expansion is conserved after the virialization, and is given by (Peebles f980; HS; White & Zaritsky f992): 

P(r)-P.i^'^ (12) 

where ri and pi are respectively the initial radius and the density, while il{z) is the density parameter at epoch z. The final 
radius, r, and the initial one, ri are connected by: 

r^Fr^ = Fri= 3-^4^ (f3) 

where Vm is the shell radius at maximum expansion, F is given in DP2000 (equation (26)) and the mean fractional density 
excess inside a given radius, 5, is: 

5 = 4/ S{y)y'dy (14) 
Jo 

In order to calculate the evolution of c = r^/a, I must calculate the inner radius a and the virial radius, Tv. The inner radius, 
a, is characterized by the condition: 

§^|. = -2 (15) 
dlog(r) 

while the virial radius, r^, is the radius within which the mean overdensity is 5v times the critical density, pc, at that redshift: 
Mv = 5v(2)pc(2)^rv(z)' (16) 
where Mv is the virial mass of the halo and the critical density pc is: 

Pciz) = po^[noii + zf + {i^no)] (17) 

The subscript indicates that the parameter is to be calculated at epoch to, Hq is the Hubble constant at to and Qo ~ ^^^,^1" ■ 
The virial overdensity, (5v, is provided by the spherical top-hat collapse model, which, for the family of flat cosmologies 
(flo + f^A = 1), gives: 

5v(z) (187r^+82y-39y^) (18) 
(Bryan & Norman 1998), where y = ^l{z) — 1 and ^l{z) is: 

niz) - ^o(^ + ^)' (19) 

In the limit z — > 00, equation (^^-([l9|) give the corresponding quantities for the SCDM model. 
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3 RESULTS 



In NFW96 and NFW97, the N-body simulations were interpreted by means of a model, a step-by-step calculation of the 
density profile that is also useful to calculate the mass and redshift dependence of the concentration parameter, c. This model 
assigns to each halo of mass Mv identified at z — zq a, collapse redshift, Zc, defined as the time at which half mass of the halo 
was first contained in progenitors more massive than some fraction / of the final mass, Mv. Lacey & Cole (1993) showed that 
a randomly chosen mass element from a halo having mass Mv, identified at redshift zq, was part of a progenitor with mass 
exceeding /Mv at the earlier redshift z with probability 



where ctq is the linear variance of the power spectrum and (5crit(z) is the density threshold for spherical collapse at redshift z. 
The collapse redshift Zc is determined setting P — 1/2. Assuming that the characteristic density of a halo is proportional to 
the density of the universe at the corresponding Zc then we have (NFW97): 



Given Mv and zq it is possible to obtain Zc from equation ( |20[ ) and Sn from equation (pl|), thus completely specifying the 
density profile. 

NFW97 model is in agreement with N-body simulations at 2 = 0, for several different cosmological models (NFW96; NFW97; 
DP2000; B99) but as shown by B99, it does not reproduce properly the redshift dependence of the halo profiles as seen in 
their simulation: it over-predicts the concentration, c, at early times, z > 1. 

In DP2000, I showed that the improved SIM model, introduced in that paper, gives good results in predicting the shape of the 
dark haloes profiles and the mass dependence of the concentration parameter, c, both in a SCDM model and in a scale-free 
universe. In that paper, I did not study the redshift dependence of the concentration parameter. Here, in order to answer this 
question, I calculated the evolution of c for two difi'erent cosmologies, namely SCDM (f2 — l,h — 0.5,n — 1; as — 0.63) and 
ACDM (r^ni = 0.3, SIa = 0.7, h = 0.7, ag = 1). The resuhs are plotted in Fig. 1-3. 

The solid line of Fig. 1 represents the expected behavior for 8 x 10^^h~^MQ haloes as predicted from the NFW97 model in a 
ACDM model. The dotted line with errorbars represents the c{z) median as obtained in N-body simulations by B99 and the 
Poisson errorbars were obtained by the quoted authors from the profile fitting procedure: after identifying a centre for the 
halo, they count particles in logarithmically spaced radial bins and assign corresponding Poisson errors based on the count 
in each bin. The short-dashed line is c{z), obtained by means of the present model. The model was normalized to match the 




(20) 




(21) 
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Figure 1. Evolution of the concentration parameter, c, in the ACDM and SCDM models of the text for a halo of 8 X 10^^ Mq. The 
solid line is the behavior of c for haloes of 8 X 10^^ Mq predicted by NFW97 analytic procedure in the case of the ACDM model of the 
text, while the dotted line with errorbars represents the evolution of c obtained in B99 simulations. The short-dashed line represents the 
prediction for c{z) by the present model. In the upper two lines (long-dashed line and dot-dashed line), I show the result of the same 
calculation for the SCDM model. Similarly to the ACDM, NFW97 model (long-dashed line) over-predicts the concentration, c, with 
respect to the present model (dot-dashed line) prediction. 

normalization of the z = relation. As can be seen from Fig. 1, the NFW97 model (solid line) overpredicts the concentration, 
c, by ~ 50% ai z — 1, with respect to the present model (short-dashed line) and by ~ 40% with respect to that of B99 (dotted 
line with errorbars). The disagreement increases with increasing z. The present model predicts coc 1/(1-1-2) in agreement with 
that of B99. In the upper two lines of Fig. 1 (long-dashed line, dot-dashed line), I show the result of the same calculation for 
the SCDM model. Similarly to the case of the ACDM model, NFW97 model (long-dashed line) overpredicts the concentration 
c with respect to the model of this paper (dot-dashed line) in which coc 1/(1 + 2:), in agreement with the analytic result and 
the simulations of B99. 

The scaling behavior of c can be explained in the same way described by B99: in the present model the scale radius, a, is 
roughly constant, then the z dependence of c = rv/a comes from the virial radius rv: aside from the z dependence of 5v, both 
in the SCDM and ACDM model rv oc 1/(1 -f z). 

Till now, I showed that the model proposed in this paper, in agreement with claims by B99, predicts a different redshift 
© RAS, MNRAS 000, §-|l4| 




8 A. Del Popolo 

dependence for the concentration c{z) with respect to NFW97. It is time to discuss the reasons behind the quoted discrepancy 
and, at the same time, the agreement of the results of the present paper with the B99 paper. 

The different result of these two models with respect to NFW97 is fundamentally due to the different way of defining the 
collapse redshift, Zc'. as previously reported in this section, NFW97 defines Zc{Mv, f, 20) as the time at which half the mass of 
the halo was first contained in progenitors more massive than some fraction / of the final mass, M^. Then the collapse redshift 
is obtained through the Press-Schechter formalism, equation (|2^) and finally the characteristic overdensity of the halo, S^, 
(which according to equation (jij) is connected to the concentration, c) is assumed to be proportional to the density of the 
universe at Zc- As a consequence of this way of defining Zc, the NFW97 prediction for c eventually goes to a constant value 
at high redshift, because Zc becomes closer and closer to zq, as can be easily found from equation (po|). In other words, at 
sufficiently high redshift the halos collapse redshift, Zc, becomes essentially indistinguishable from the redshift, zo, at which 
they are analyzed and their concentrations tend to a constant. 

Things go differently with respect to the previous discussion, both in the present paper and in that of B99. In the present 
paper, I use the SIM, which, as it is well known, allows one to establish, for a given power spectrum, a relationship between 
the mass Af and its formation (collapse) epoch, Zc (see Peebles 1980, Section 19; Gunn & Gott 1972; Avila-Reese, Firmani & 
Hernandez 1998). In fact, given a density perturbation, a shell with initial comoving radius ri and mass M = ^pi,rf{l + Si), 
will expand to a maximum radius rm, given by equation ( |l^ ) or by rm — ri/Si in an Einstein-de Sitter universe, and then 
collapse will occur. The time of maximum expansion in an Einstein-de Sitter universe is given by: 

tc TT 1 

where tc is the collapse time and Si is the mean fractional density excess measured at time ti. In other words, equation ( |l^ ) and 
equation (^) tell us that the collapse redshift z^ is a function only of the virial radius (and then of the virial mass, Mv), and 
is independent of zq, in agreement with the assumption made by B99 in their equation (8). As a consequence, the condition 
Zc — Zo, encountered in the NFW97 model at high redshift, implying that c[z) tends to a constant value, is no more present in 
the model of this paper. This implies that the slope of the relation c-z remains constant. In conclusion, c{z) oc 1/(1 + z) and c 
never tends to a constant. The same result is found in B99 model. I also would like to remark that in the present model and in 
that of B99, it is possible, at high redshift, that Zc < 20 if a halo is much more massive than the characteristic mass, Af, Q at 
that epoch, that is, the characteristic collapse time for the subclumps may take place in the future. This strange result means 
that the collapse of the sublumps happened only because the large scale structure, within which they were sitting, collapsed, 

* I recall that as shown in DP2000 for the adopted normalization of the CDM spectrum Af* = 3 X lO^^AfQ 
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Figure 2. Evolution of the concentration parameter, c, in the SCDM model of the text as function of the mass and for different redshifts. 
The solid lines represent the prediction of the model of this paper at 2: = 0, z = 2, 2: = 4. The dotted lines represent the prediction of 
NFW97 for the same redshifts and for the values of the parameters, constant in time, / = 0.01 and C = 3 X 10*, defined in NFW96 
and NFW97, while the dashed lines represent the predictions of B99 model for the values of the two parameters, constant as a function 
of both z and mass, F = 0.01, K = 3.8, defined in their paper. The filled hexagons represent the concentration c at 2 = obtained by 
NFW96 and NFW97 N-body simulations. 



or, in other words, the collapse of the subclumps happened at the same time as the collapse of the whole halo itself. In this 
case the halo has a very low concentration. This implies that haloes much larger than Af, will always be of low concentration, 
because their subclumps did not collapse much earlier than the halo itself. We may also add that the concentration declines 
dramatically with z because haloes at high redshift are just beginning to form and are therefore more diffuse. 
A possible question that may arise at this point is why the NFW97 model works well at z = and not at higher redshift. 
In fact at z — 0, the NFW97 model correctly predicts the mass-density relation obtained from N-body simulations and the 
result is also in agreement with B99 and DP2000. The answer to the previous question is that the problem of NFW97, seen 
at z > 0, is no more present at z = because at this redshift, the extended Press- Schechter formula, that NFW97 used to 
determine Zc, never gives 2:0 — Zc, so the problem of the wrong prediction of c(z), seen at high redshift, is no more present 
at z = 0. In order to further clarify the previous discussion between the differences in the model of this paper and NFW97, 
I calculated the concentration c as a function of mass, Af, in a CDM model for different values of redshift in the case of 
© RAS, MNRAS 000, §-|l4| 
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NFW97 model, in that of B99 and in the model of this paper. The result is shown in Fig. 2. The solid lines represent the 
prediction of the model of this paper at three different redshifts, z = 0, z = 2, z = 4. The dotted lines represent the prediction 
of NFW97 for the quoted redshifts, while the dashed lines are the predictions of B99 model. The filled hexagons represent the 
concentration c at z — obtained by NFW96 and NFW97 in their N-body simulations. Fig. 2 shows that for a fixed mass and 
at z = 0, the evolution of c as a function of mass AI is, as known, succesfuUy predicted by NFW97 model for the values of the 
parameters, constant in time, / = 0.01 and C — 3x lO**, defined in NFW96 and NFW97. B99 model (for the two parameters, 
constant as function of z and mass, F — 0.01, K = 3.8, see their paper for a definition) and the model of this paper reproduce 
the z = results of NFW97 quite weU over the range jf- ~ 0.01 — 100 (being as previously reported Af, = 3 x 10^^ Mq). 
A direct comparison of the three models however shows that the NFW97 model prediction is shallower than the other two. 
The same situation was present in the ACDM model, as shown in Fig. 2 of B99. Also in this case, NFW97 predicted slope for 
c is shallower than B99 model which reproduce better the simulations data. At 2: = 2 the discrepancy between the NFW97 
and the other two models is evident. Fig. 2 evidently shows that, going from z = to z = 2, the slope of the relation c — M, 
predicted by B99 and the model of this paper, remains roughly constant and moreover gives lower values of c with respect to 
the NFW97 model. The figure also shows that the NFW97 model for 2 = 2 have a smaller slope with respect to the z = 
case. The situation now described is even more evident in the 2 = 4 case. This situation was expected because, as previously 
quoted in the precedent discussion, with increasing 2, 2c —> 20, and one expects that c goes to a constant, in fact we see that 
the slope of the c relation predicted by the NFW97 models decreases going from 2; = to 2 = 4. In other terms, the B99 model 
and the model of this paper predict that c{z) oc 1/(1 + 2) and that the relation never tends to a constant for a particular 
high redshift, which is the case for the NFW97 model. As previously quoted, the reason why NFW97 model gives a succesful 
description of the c — M relation in the case 2: = is due to the fact that in this case never happens that Zc — zq. 
Finally, I want to remark that although B99 model gives predictions for the c(M, 2) relation in good agreement with N-body 
simulations and with the model of this paper, there is a fundamental difference between these two models. B99 model makes 
the 'ad- hoc' assumption that Zc depends only on mass and not on 20, in order to explain the results of the N-body simulations 
performed in the same paper. The model of the present paper is a direct consequence of SIM, without additive assumptions. 
Now I want to discuss some consequences of the results obtained for c(z) on galaxy modelling at high redshift. The discussion 
is intended to gain only a qualitative understanding of how c{z) may affect some features of galactic formation (a more 
complete discusion is given in B99). 

As pointed out by B99, one expects that the described behavior of c{z) should have a large impact on galaxy modelling at 
high redshift and for interpreting high redshift data (e.g., evolution of the TuUy-Fisher relation (Vogt et al. (1997) and the 
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Figure 3. Disk size, Rd, as a function of z for the SCDM and ACDM model introduced in tlie text assuming that Mv = 8 X 10" Mq . The 
solid and dotted line represents for the SCDM model, obtained from the B99 fitting formula (their Eq. (7)), respectively assuming 
that the evolution of c follows NFW97 model (solid line) and the model of this paper (dotted line). The long-dashed and short-dashed 
line represents for a ACDM model obtained from the B99 fitting formula respectively assuming that the evolution of c follows NFW97 
model (long-dashed line) and the model of this paper (short-dashed line). 

nature of Lyman Break Galaxies (Steidel et al. 1996)). According to the standard picture of galaxy formation, structures grow 
hierarchically from small, initially Gaussian density fluctuations. Collapsed, virialized dark matter haloes condense out of the 
initial fluctuation field. Gas associated with such dark haloes cools and condenses within them, eventually forming galaxies. 
In this scenario, the growth of the dark matter haloes is not much affected by the baryonic components, but determines how 
they are assembled into nonlinear units. The halo density profile determines many of the properties of galaxy disks, e.g. their 
size and surface brightness. In order to show the effect of c evolution on disk properties, it is possible to use a fitting formula 
given by Mo, Mao & White (1998) based on several assumptions about the halo and disc make-up and depending on some 
free parameters, or more easily the fitting formula introduced in B99 (Eq. (7)): 

= 5.r/.-Kpc (^^5^^) [1 + (.(.)/3.73)° -] - (23) 

which for 1 < c < 50 is good within 1%. In Fig. 3, I plot the z dependence of the disk size, Rd, in the SCDM and ACDM 
model, assuming that Mv = 8 x 10^^ Mq, as assumed also in Fig. 1. In both cases the lower concentration of haloes at high 
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12 A. Del Popolo 

redshift produces a disk size larger than that obtained using the NFW97 model for the evolution of c. In fact in the ACDM 
model taking account of the c — z dependence found in this paper, I find the short-dashed line, while the long-dashed line is 
obtained with the NFW97 c — z dependence. A similar situation is found in the SCDM model: using the c — z dependence of 
this paper, I obtained the dotted line while using that of NFW97 I obtained the solid line. I want to stress that the last result 
must be taken with caution, since it strongly dependes on the precision of the fitting formula used and on its free parameters 
(if it contains any). In fact comparing Fig. 1 and Fig. 3, it is immediately evident that, while the concentration found in the 
model of the present paper differs from that predicted by NFW97 evolution quite a lot, the differences is reduced when we 
deals with _Rd. If I had used the fitting formula of Mo, Mao & White (1998), to predict the variation of _Rd with redshift, I 
should have found that small variations in the parameters of that model, at a level of 10%, would have produced change in Rd 
larger than that due to the redshift evolution of c. So the validity of the last result, and that regarding the surface brightness, 
I{z), described in the following, is entirely based upon the reliability of the B99 fitting formula. 
Since the surface brightness, I{z), scales as: 

I{z) cc Rf{z) (24) 

both in the SCDM and ACDM model, we may expect that this parameter is also modified by the strong evolution of c 
with redshift. Using Fig. 3 and equation ( p^ it is easy to find that according to the model of the present paper disks are 
dimmer with respect to NFW prediction and that the effect increases with increasing z. Besides the two examples that I have 
just given, the strong evolution of the concentration parameter, c, affects several other properties of disks and, as previously 
reported, has some implications for galaxy modelling. For example, the shape of a disk rotation curve depends, among other, 
on the concentration, c, of its halo. Since more strongly peaked curves are found in more concentrated haloes and evolution 
produces a reduction of the concentration of haloes one expects that evolution tends to produce less peaked rotation curves. 
Further implications of the c{z) dependence shown in this paper are discussed in B99. However, as stressed by the quoted 
author, much work remains to be done in order to have a deeper understanding of the implications of the results obtained in 
this and their paper on galaxy modelling. 



4 CONCLUSIONS 

In this paper I studied the evolution of the concentration parameter, c{z), for fixed mass haloes, in the SCDM and ACDM 
model by means of the improved SIM introduced in DP2000. The results of the paper can be summarised as follows: 
a) Both in the SCDM and ACDM model the evolution of c{z) is much stronger than that expected from the NFW97 model: in 
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the model of the present paper c{z) cc 1/(1 + 2), while NFW97 model overpredicts the concentration c by 50% ai z — 1, with 
respect to the model of the present paper, and by 40%, with respect to that of B99. NWF97 predicts that the concentration 
c tends to a constant value for large z, so the quoted disagreement increases with increasing z. 

b) A comparison of the results of the present paper with the high-resolution N-body simulations of B99 shows a good agree- 
ment both for the SCDM and ACDM model. 

c) The different redshift dependence of c{z), obtained in the present paper and B99, with respect to NFW97, is fundamen- 
tally due to the different way of defining the collapse redshift, as a consequence of the way of defining Zc in NFW97, at 
sufficiently high redshift the halos collapse redshift, Zc, becomes essentially indistinguishable from the redshift, zo, at which 
the halos itself are analyzed and their concentrations tend to a constant. 

d) The present model has an advantage with respect to B99, namely: B99 results are obtained assuming the 'ad-hoc' as- 
sumption that Zc depends only on mass, and not on Zo, while the results of the present paper are direct consequence of SIM, 
without additive assumptions. 

e) The result has important consequences on galaxy modelling at high redshift and on some disks characteristics. In particular 
the disk size obtained in the present model is larger than that obtained using the B99 fitting formula together with NFW 
prediction for c{z). The reverse is true for the disk brightness. 
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